A NOTE ON THE DYNAMICS AND 
THERMODYNAMICS OF DISLOCATED CRYSTALS. 

RUSLAN SHARIPOV 

Abstract. The dynamics and thermodynamics of dislocated crystals are studied 
within the framework of the nonlinear theory of elastic and plastic deformations. 



1. Introduction. 

In the series of papers [1 4] the kinematics of crystalline media in the presence 
of dislocations was studied. Here we consider the dynamics and thermodynamics 
of a crystalline medium. This paper continues the series of papers [1-4]. Therefore, 
we do not provide a large introductory section in it. For more details the reader is 
referred to the previous papers. 

2. Specific free energy function. 

By analogy to [5], we shall construct the theory on the base of the specific 
free energy function / (the amount of free energy per unit mass of a medium). 
The deformation state of a crystalline medium in the presence of dislocations is 
described by the incompatible distorsion tensor T. This is a double space tensor 
with one upper index 1 associated with the Burgers space B and with one lower 
index associated with the real space E. The other physically important parameter 
is the tensor of the Burgers vector density p. It is also a double space tensor with 
two indices. Like in the case of T, the upper index of p is associated with the 
Burgers space, while its lower index is associated with the real space. The easiest 
way of choosing the specific free energy function / is to write it as 

/ = /(T,f,p), (2.1) 

where T is the temperature. However, we shall not write / in this simplest way (2.1). 
The matter is that T and p contain some excessive amount of purely geometric 
information absolutely inessential for describing the thermodynamics of a medium. 
We need to make some efforts to get rid of this inessential part of information in 
the tensors T and p. 

The Burgers space B, as it was defined in [1], is a copy of the real space E. It is 
assumed to be filled with the ideal (dislocation-free and non-distorted) crystalline 



It is implicitly assumed that some coordinate systems in B and E are chosen and T is repre- 



sented by its components T? respective to those coordinate systems. 
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medium (see Fig. 2.2) identical to that we consider in the real space. The tensor 
field T determines a linear mapping at each point p of the real space: 

f : T p E^M. (2.2) 

Here T p K is the set of vectors attached to the point p € E. In geometry, when E 
is treated as a manifold, such a space is called a tangent space. Both spaces E and 
B are equipped with the Euclidean metrics g and g respectively. The difference is 
that the Burgers space B in (2.2) is treated as a linear space, while E is treated as 




Fig. 2.1 Fig. 2.2 

a point space (or a manifold). The mapping (2.2) is non-degenerate: detT ^ 0. 
Hence, one can consider the inverse mapping 

S: B T p E. (2.3) 

The inverse mapping (2.3) is given by the inverse incompatible distorsion tensor S. 
This is also a double space tensor with two indices: its upper index is associated 
with E and its lower index is associated with B. 

Suppose that the point p is fixed. The rotation of the crystalline body as a 
whole about the point pgE does not change its physical state (see Fig. 2.1). Such 
a rotation is expressed by the following transformations at the point p: 

3 3 
a— 1 a— 1 

Here 0£ are the components of some rotation matrix 1 . Since the physical state 
of a crystalline body is invariant under the transformations (2.4), so should be the 
free energy function (2.1). Let's define the following tensors G and R: 

3 3 3 

^ = E E 9c.fi S? S?, R] = Y,P« ( 2 - 5 ) 

a=l [3=1 a=l 



1 Note that the transformations (2.4) are different from those in section 2 of [3] and from those 
in formulas (3.19), (3.20), (3.21) in [4]. 
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Both indices i and j in (2.5) are associated with the Burgers space B. Nevertheless, 
G and R are double space tensors because their components are functions of the 
coordinates of a point peE and of the time variable t: 

G l3 =G l3 (t,y\y 2 ,y 3 ), Rj = Rj(t,y\y 2 ,y 3 ). (2.6) 

Note that the tensors G and R are different from the purely real space tensors G 
and R introduced in [1] and in [3]. Their components are given by the formulas 

3 3 3 
i=l j=l i=l 

The tensor G given by the first formula (2.7) is called the elastic deformation 
tensor. As for R, it is the purely real space representation of the Burgers vector 
density tensor p. 

Now let's return back to the tensors G and R in (2.5) and find that they are 
invariant with respect to the transformations (2.4). Therefore, we can use G and 
R as independent variables instead of T and p in (2.1): 

/ = F(T,G,g,Rp). (2.8) 

Note that the function F has no explicit dependence on the coordinates of a point 
p £ E. This function represents the refined properties of a crystal in its homoge- 
neous state (see Fig. 2.2). It should be tabulated on the base of experiments. Note 
also that, apart from G and R, we introduced the additional arguments g and p in 
(2.8), where p is the density of a crystal and g is the metric tensor of the Euclidean 
metric in the Burgers space B. The constant tensor g does not change the time 
and coordinate dependence of / determined by (2.6) and by the scalar functions 

T = T(t,y 1 ,y 2 ,y 3 ), p = p(t,y\y 2 ,y 3 ). 

Having written (2.8) in place of (2.1), we have reached one of our goals — we 
have got the rotationally invariant free energy function. The experimental data 
should go to the theory through the function F. However, for the sake of beauty, 
in the formulas below we need to express / through some purely real space tensor 
fields. The tensor R is expressed through R as follows: 

3 3 
as=l 0=1 

The tensor G cannot be expressed through G, but it is expressed through the 
metric tensor g of the real space E according to the first formula (2.5): 

3 3 
a=l /3=1 
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Similarly, the metric tensor g in (2.8) is expressed through G as follows: 

3 3 
o=l 0=1 

This formula (2.11) is inverse to the first formula (2.7). Substituting (2.9), (2.10), 
and (2.11) into (2.8), we can write (2.8) in the following form: 

f = f(t,p,T,g,G,TL,p). (2.12) 

We do not keep S and T in (2.12) as independent variables since they are not 
purely real space tensor fields. Their presence is reflected by the arguments t and 
p, where t is the time variable and p is a point of the real space E. In a coordinate 
presentation of / its argument p is replaced by the variables y , y 2 , y 3 . 

3. Specific free energy function 
as an extended scalar field. 

The specific free energy function (2.12) is a scalar function with three tensorial 
arguments g, G, and R. In this form it fits the definition of an extended scalar 
field (see definition 4.1 in [6]). The theory of extended tensor fields was especially 
derived in [6] for to use it in the present paper. The extended scalar field / in 
(2.12) is not an arbitrary extended scalar field. It is produced from the function 
F in (2.8). For this reason it satisfies some partial differential equations written in 
terms of covariant and multivariate derivatives defined in [6]. 

Theorem 3.1. The specific free energy function f of a dislocated crystalline medi- 
um is an extended scalar field satisfying the differential equation 

3 3 3 

V 7 /-EEE 9mp + g am ) - 

m=l a—1 3— 1 
3 3 3 

- E E E G m p + G am ) + (3.1) 

rn — 1 a—1 (3—1 

+ E E E E v £ w {^m k - z™ Rt) , 

n—1 m—1 a—1 (3—1 

where Z^- are the components of the tensor field Z introduced in paper [3] : 

3 3 / Qrjvtn 3 \ 

^ = E^ V «T™ = E*M^-E r ^ 1 ■ (3 . 2) 

m=l m=l \ ^ )j = l / 

Proof. First of all let's calculate the multivariate derivatives V a/3 [2]/, V a/3 [3]/, 
and Vp[4]f defined due to three tensorial arguments in (2.12): 

V*[2]f=V=£±-]*S?§!, (3.3) 
a 9°P i=i j=1 ddj 
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Then we calculate the covariant derivatives of the extended scalar field (2.12) ap- 
plying the formula (13.15) from [6]: 

v-rf - 1£ + E E E ^ + r ^ + 

+ EEE ( r ;« ^ + r ^ + (3.6) 

m=la=l/3=l C«a/3 
3 3 3 y> 

+ E/ E^ EZ (^7" ^™ ~ ^ m ) anP ' 
m=la=lp=l OH a 

Writing (3.6), we assume that y 1 , y 2 , y 3 are some curvilinear coordinates in the 
real space E and T^- are the components of the metric connection for the Euclidean 
metric g in E. Now let's calculate the partial derivatives df jdy 1 : 

df 3333 , a dSf\ 

1=1 j=l a =l 0=1 \ y y / 

The matrices S and T are inverse to each other. Applying this fact and the equalities 
(3.3), (3.4), (3.5) to the above formula, we obtain 

df f f f f df (~ m dfs ,™ 5f >" 

^--l^l^l^Pd^-A S n d^9mP + S n —g, 
" n=\m=\a=\0=\ ^ ap \ y " 



E\ ^ \ " \ " ®f ( cm ®T£ A , cm ® T & A 



+EEEE^PnT M <-sTTv< ■ 

Now, substituting this expression for df jdy 1 into (3.6) and taking into account 
(3.2), we derive the required equality (3.1). The theorem is proved. □ 
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Note that in (3.2) and (3.6) are the components of the symmetric Euclidean 
connection T associated with the standard Euclidean metric g in the real space E. 
These quantities are equal to zero in Cartesian coordinates and arise only if one 
chooses curvilinear coordinates y , y 2 , y 3 . Let's remember that there is another 
connection T in the theory. Its components are given by the formula (2.3) in [4]: 

f*=r£. + £*. (3.7) 

This connection (3.7) is associated with the tensor field G treated as a metric. If 
we replace T by T in (3.6), then the equation (3.1) reduces to the following one: 

%f = 0. (3.8) 

Remark. Note that both V 7 and V 7 in (3.1) and (3.8) are not standard covariant 
derivatives. They are so-called spacial covariant derivatives introduced in [6]. 

Now let's calculate the time derivative df /dt considering / as an extended scalar 
field (2.12) and thus taking T, g, G, R, and p for independent variables: 



df _ **** df (~df S ,<w^" \ 

dt - 2^2^ 2^2^ dg I b n qI 9m0 + b n g am I 

n=l m=l a=l f3=l a ' \ J 

E\ ^ \ \ " &f ( dm ~ ~ rn dTp ~ \ 

I °n -qT + O n -77- Lr am I + 

n=l m=l a=l 0=1 a(j afi \ ) 

n=lm=l a=l (3=1 UJX a \ / 

The components of the time derivative dT/dt are given by the formula 

Q r J^ n 3 3 

-gf = -J f 3 E V ^ ^ ~ E u ' (3-9) 

i=l i=l 

The formula (3.9) is derived from the formulas (4.4) and (4.5) in [2]. These two 
formulas are valid only under the assumption that the conjecture 4.1 in [2] is valid. 

Remark. Here we do not discuss the conjecture 4.1 from [2]. However, one should 
note that all results below in this paper are obtained under the assumption that 
this conjecture is valid. 

The quantities jp in (3.9) represent the double space tensor j , it is called the 
tensor of the Burgers vector flow density. The quantities v l are the components 
of the velocity vector v describing the substance flow in the medium. Multiplying 
(3.9) by S™ and summing over the index n, we derive the following equality: 

3 ST n 3 
E S™ = -jp - Vpv m - E v'ZR. (3.10) 

n— 1 i—1 

Here Jq 1 are the components of the tensor J. This tensor is the purely real space 
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version of the tensor j . It is also called the tensor of the Burgers vector flow density. 
The components of the tensor J are given by the formula (2.9) in [3]: 

3 

= (3.11) 

n=l 

Apart from (3.11), let's recall also the formula (2.10) from [3]. It can be written as 

3 

0™ = + £ yl ^"h - 33i)- ( 3 - 12 ) 
»=i 

The tensor 8 with the components (3.12) is called the tensor of the rate of plastic 
relaxation. Applying (3.12) to (3.10), we derive the following formula: 

3 QT n ^ 

£ s™ -JL = v p v m - £V zfr (3-i3) 

n—1 z—1 

Now let's apply (3.13) in order to transform the above expression for df /dy 1 : 

f-* p 3 3 3 3 

i+E l '' V '/ = -EEE VQ/3 [ 2 ]/ K 9mp - V a v m g mp + 

z—1 rn — 1 a — 1 j3—l 

3 3 3 

+ ^" g am - Wpv m g am ) - £ £ £ V^[3]/ (C G m p - 

m=la=l/3=l (3.14) 
3 3 3 

- V qU " 1 G m/3 + 0^ G Qm - V^™ G am ) + E E E V W X 

m =l a =l /3=1 

The result obtained in the formula (3.14) is formulated as the following theorem. 

Theorem 3.2. The specific free energy function f of a dislocated crystalline medi- 
um is an extended scalar field satisfying the differential equation (3.14). 

4. Balance equations. The dynamics and thermodynamics. 

Like in [5], here we shall develop the theory on the base of the balance equations. 
The first balance equation is very simple. This is the mass balance equation: 



d tP + J2^k(pv k )=0. (4.1) 

fc=i 

The momentum balance equation looks not much more complicated than (4.1): 

3 

a t ( P ?/) + £ v ar fe = /\ (4.2) 



k=l 
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However, the whole complexity of this equation is hidden in the tensor II. As for 
the quantities p in the right hand side of (4.2), they are the components of the 
vector f . They define the density of volume forces acting upon the medium. The 
components of the tensor II arc given by the following formula: 

IT fe = pv l v k - a ik - a ik . (4.3) 

Here a %k and u lk are the components of two stress tensors er and er respectively. 
The first of them er is the regular stress tensor, we shall study it in more details 
a little bit later below. The second tensor er in the formula (4.3) is the viscosity 
stress tensor. Its components are given by the formula 

^ fc = ^EE^( V ^ + V ^), (4-4) 

j=l q=l 

where rj i q are the components of the viscosity tensor r\. 

The third balance equation expresses the energy balance in a medium. This 
equation is written in the following form: 

^ ' k=i 

Here e is the specific inner thermal energy. We shall discuss it a little bit later 
below along with the regular stress tensor er. The quantities w k in (4.5) are the 
components of the vector w. They are given by the formula 

w k = p\^_ v k +p£V k_j2 Vi a ik - ~^ - E v * r * lk - ( 4 - 6 ) 

i—1 i—1 i—1 

The vector w with the components (4.6) is responsible for the energy transport in 
a medium. It is called the energy flow density vector. The first two terms in the 
right hand side of (4.6) correspond to the energy transported with the mass flow. 
The next two terms represent the work performed by the stress forces. And the 
last term describes the heat conductivity phenomenon. The quantities x lk are the 
components of the heat conductivity tensor. 

The quantity e in the right hand side of (4.5) describes the energy sources within 
the bulk of a medium. In the simplest case the energy is produced by the work of 
the volume force f . In this case we can write 

e = E u */ i > ( 4J ) 
i=i 

where /* are the same as in (4.2). They are the components of the vector f . 

Our further goal is to specify the above three balance equations (4.1), (4.2), 
and (4.5) for the case of a crystalline medium with dislocations. The specific inner 
thermal energy e is related to the specific free energy / as follows: 



e = Ts + f. 



(4.8) 
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Here s is the specific entropy. It is derived from / by means of the formula 

In this form the formula (4.9) is applicable to (2.8) and to (2.12) as well. However, 
applying it to (2.12), we can use the notations introduced in [6]: 

a = -V[l]/ = s(t,p,T,g,G,R,p). (4.10) 
Applying (4.9) to the function (2.8), we obtain the following equality: 

s = -^ = S(T,G,g,R,p). (4.11) 

The functions S and s in the right hand sides of (4.11) and (4.10) are related to 
each other in the same way as the functions F and / in (2.8) and (2.12). For this 
reason we can write the differential equations 

3 3 3 
m=l a =l /3=1 

-EEE VaP W s d ^ + ^ G am ) + (4.12) 

m=l a=l [3=1 
3 3 3 3 

+/Z/Z/ZY, V W k - Ri) , 

n— 1 m— 1 a— 1 /3— 1 
3 3 3 3 

E" <ViS = -EEE VQ/3 [ 2 i s ( ™ ^ - v ^ m ^ + 

2—1 m— 1 a— 1 /3— 1 

3 3 3 

+ 0^ g am - V p v m g am ) - E E ^t 3 ] 5 (C 6 '"/3 " 

m=l a=l /3=1 

3 3 3 

V Q f" 1 G ro /3 + Gam ~ VpV™ G a m) + E E E V £ M SX 

m =l a= l /3 = i 



ds 

dt 



(4.13) 



x (C i?" - V m / iC - C < + V a v m R 

Theorem 4.1. The specific entropy of a dislocated crystalline medium s is an 
extended scalar field satisfying the differential equations (4.12) and (4.13). 

It is known that S is a monotonic function of T in (4.11). The same is true for 
s in (4.10). Moreover, we have the thcrmodynamical inequalities (sec [7]): 

dS ds , 

W > 0, W > 0. (4.14) 

This means that wc can introduce the inverse functions 



T = T(s,G,g,TL,p), 



T = T(t,p,s,g,G,TL,p). (4.15) 
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Applying the inequalities (4.14) to the functions (4.10) and (4.11), we derive 

c = r|i>o. 

Here c is the specific heat capacity of a medium for the case where the density p 
and the other arguments t, p, g, G, R, G, g, R in (4.10) and (4.11) are constants. 

Let's substitute the temperature expressed in the form of the function (4.15) 
into (4.8). As a result we get two functions 

£ = E(s,G,g,&,p), £ = e(t,p,s,g,G,n,p). (4.16) 

Passing from (2.8) and (2.12) to the functions (4.16), we perform the so-called 
Legendre transformations. The functions (4.15) are related to each other through 
(2.9), (2.10), and (2.11). The same is true for the functions (4.16). The second 
function in (4.16) is an extended scalar field, though with the slightly different set 
of arguments than in (2.12). From (2.9), (2.10), (2.11), and (4.16) we derive 

^me=^ = ±±^S?Se, (4.17) 

V^[3] £ =4^ = X^-^£f5f, (4.18) 
dG a p 2—i j— i dgij 

V|[4] £ =^=^^^;^^. (4.19) 
Applying (4.17), (4.18), (4.19), and the formula (13.15) from [6], we get 

3 3 3 

V 7 £ = - E E E ^^t 2 ] 5 fa, 9mp + Z™ p g am ) 

m=l a=l /3=1 
3 3 3 

- E E E ^ + Z™ G am ) + (4.20) 

m=l a=l j3=\ 

+ E E E E v l W £ {^m k - z™ a Ri) , 

n— 1 m— 1 a— 1 /?— 1 

Like the equation (3.1), the equation (4.20) can be simplified by means of (3.7): 

V 7 e = 0. (4.21) 

Note that the remark following the formula (3.8) is valid for the covariant derivatives 
V 7 and V 7 in the above formulas (4.20) and (4.21) as well. 

Theorem 4.2. The specific inner thermal energy s of a dislocated crystalline medi- 
um is an extended scalar field satisfying the differential equation (4.20), or the 
differential equation (4.21) equivalent to it. 
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From (4.16), (4.17), (4.18), (4.19), and from (4.20), applying (3.13), we obtain 

3 3 3 3 



i=l m=la=l/3=l 

3 3 3 

+ 9am ~ VpV™ 9am ) " E E E Va0 W £ - 

m=lo=l/3=l (4.22) 
3 3 3 

- V a V m G mf} + 0^ G am ~ Vf}V m G am ) + E E E V /3 Q M £X 

m=l a=l /3=1 

x K K - v ro ^ K - C i& + v Q « m O . 

Theorem 4.3. The specific inner thermal energy e of a dislocated crystalline 
medium is an extended scalar field satisfying the differential equation (4.22). 

The equalities (4.17), (4.18), and (4.19) should be completed with the following 
one being analogous to the equality (4.9): 

T=^. (4.23) 
us 

The equality (4.23) applies to both functions (4.16) and to both functions (4.15): 

<9e(s,G,g,R, p) 



T(s,G,g,R, p) 



d 



s 



rru r< v) ^ de(f,p, s, g, G, R, p) 

T{t,p,s,g, G,R,p) = — = V[lJe. (4.24) 

Now let's return to the balance equation (4.5). Substituting (4.6) and (4.7) into 
(4.5), then taking into account other two balance equations (4.1), (4.2) and the 
formula (4.3), for the function e we derive 

3 3 3 

pd t s + E ^ 9fc£ = E E ( VfcW * + + VjfeCWTx*)) . (4.25) 

fe=l i=l k=l 

(see more details in [5]). All tensor fields in the balance equations (4.1), (4.2), 
(4.5) and in the formulas (4.3), (4.4), (4.6), (4.7), (4.25) are treated as regular 
(not extended) tensor fields and covariant derivatives Vfe and V, are understood as 
regular covariant derivatives different from those in left hand sides of (4.20) and 
(4.22). The partial derivative dt£ is expressed through de/dt in (4.22) as follows: 

Bf 3 3 

dte = — + V[l]e d t s + E E VaP W e d ^ + 

" = 1/3=1 

4.26) 

3 3 3 3 K ' 

+ E E vaP ^ £ d ^ + E E W]e dtRi + V[5]e d tP . 

a=\f3=\ 0=1/3=1 
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Similarly, for the partial derivative d k e in (4.25) we can derive the formula express- 
ing it through the covariant derivative (4.20): 

3 3 

d k e = V k e + V[l]e d k s + E E V a/3 [2] £ V fe g af} + 

a=l 0=1 

(4.27) 

3 3 3 3 y ' 

+ E E V ° P W £ v ^ + E E v » v ^ + d *p- 

a=10=l a=l 0=1 

In deriving (4.27) we used the formula (17.9) from [6]. Now let's use (4.24), (4.26), 
and (4.27) in order to calculate the following expression: 

d t e + E v k d k e = ( ^ + E v k V k s J + T (d t s + ^ s J + 

fe=i V fc=i / V fc=i / 

3 3 / 3 \ 3 3 

+ E E val3 W £ d ^ + E vk ^G aP + E E v /» x ( 4 - 28 ) 

0=10=1 V k=l ) at=\0=l 

x ^flg + E «* V ^ + V[5] £ ^> t p + E 9 k pj . 

The terms containing dtg a p and V k g a p do vanish due to the following equalities: 
d t g a p = 0, V fc5Q(3 = 0. (4.29) 

The first equality (4.29) is obvious: the metric tensor g is a geometric equipment of 
the real space E, its components are constants. The second equality (4.29) expresses 
the concordance condition for the metric and connection (see [8] and [9]). 

Let's apply the formula (4.22) and the first balance equation (4.1) in order to 
transform the equality (4.28). As a result we write it as follows: 

3 / 3 \ 3 3 

d t £ + E^ £ = T U s + ^ fe d fcS +EE VQ/3 t 2 ] £ x 

k=l \ k=l ) a=X0=l 



(3 \ 3 3 

V a v p + V /3Wq -E K 9mp + 9? g am ) + E E VQ/3 t 3 ] £ 
m=l ) a=l 0=1 

/ 3 3 3 

X 

\ k=l m—1 m—1 



d t G aP + E "* VkG a ,3 + E V « v " 1 + E V /3 v " 1 G ™ 

V, A;— 1 m—1 m—1 

3 3 \ 3 3 / 

- E c G ^ E °p n G ^ + E E v ^w e ^ -+ 

m=l m=l / o=l /3=1 \ 



(4.30) 



E « fc v fc i?f - e Vm^ ic + £ ^ m R ™ + E C R ™ 

k—1 m—1 m—1 m—1 

3 \ 3 



v [ujt vfe'<' 

fc=l 
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Now it is worth to remember the differential equation for the elastic deformation 
tensor G. It was suggested empirically for plastic media in [5]. Later it was derived 
for dislocated crystalline media in [1] (see formula (3.8) over there): 



dG 



dt 

r=l ^ r=l ^ r=l ^ ^ 

+ E ®k + E Oq- 

r—1 r— 1 

Comparing (4.31) and (4.30), we find that the term with V Q/3 [3]e in (4.30) vanishes 
due to the equation (4.31). As a result we get 

3 / 3 \ 3 3 

dte + E v k d k e = T ld t s + E v k d k s + E E x 

fe=l V fe=l / a=l 0=1 

(3 \ 3 3 

Vcvp + Vpv a ~Yl (C 9mP + O 7 n g am ) + E E V ^ M £ x 
m=l ) a=l 0=1 , . 

( 3 3 3 \ • I 

mi + E yk v * fl a - E V ^ vf3 R ™ + E v ^ m R ™ + 
k—1 m—1 rn—1 

3 3 \ 3 3 

+ E °™ k - E e » < + E E ^ p 3 Q/3 

m— 1 m — 1 / a— 1 /3— 1 

In the next step we should recall the differential equation that determines the time 
evolution of the tensor R (see (4.11) in [3] or (4.1) in [4]): 

f)R0 3 3 3 3 

= E J ™ R « + E c + EE" p *s? - 

m—1 m—1 m—1 p—1 

4.33 

3 3 3 3 3 3 3 v ' 

- E E E ^™ ^ mrs Vr jf - E E E E 3— ^ mrs ^ J * ■ 

m=l r—1 s— 1 m—1 r—1 s — 1 p—1 

Let's apply the relationship (3.12) to the term Jj^ in the first sum of (4.33): 

QD/3 3 3 3 3 

= - E < R ™ + E c + EE vP k - 

m—1 m—1 m—1 p—1 

3 3 3 3 3 3 3 ( 4 ' 34 ) 

-EEE3™ V ^ E E E E 9cm u>™ Z% Jf, 

m—1 r—1 s— 1 m—1 r—1 s— 1 p—1 

Then let's recall that the formula (3.12) can be written in a different form: 

3 3 3 

J," = -OF + E E E V g™ RP (4.35) 

r— 1 7—1 n— 1 
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(see formula (2.13) in [3]). Before substituting (4.35) into the equation (4.34) let's 
recall the following purely algebraic identity (it was already used in [2] and [3]): 



Y ^ s u mrs = 8™ S; - <5™ 5 r T . (4.36) 

8=1 

Here <5™, 5^, 5™, and 6^ are Kronecker <5-symbols. Now substituting (4.35) into the 
last term of (4.34) and taking into account (4.36), we derive 

EEEEs™ - mrs % J * P = -EEEE *«» - mrs ^ P 0* + 

m— 1 r— 1 5 — 1 p— 1 m— 1 r—1 s—1 p—1 

+ EEEEEEE ^ " mrs % ^ * 9 ™ k = 

m— 1 r— 1 s—1 p—1 r— 1 7— 1 n— 1 

3333 333333 

= - E E E E ^ - mrs ^ P + EEEEEE ^ >< 

m— 1 r— 1 s—1 p—1 m — 1 r— 1 p—1 r— 1 7— 1 n — 1 

3 3 3 3 

x (C 5; - £™ 5 r T ) Z% v< g™ RP = Y, E E E 9am u mrs Z% X 

rn — 1 r— 1 s—1 p—1 

3 3 3 3 3 

x«r+EE" r ^ P ^ - E E E w « ^ s rn ^ ■ 

r—1 p—1 r—1 p—1 n— 1 

Applying this result to the equation (4.34), we can write this equation as follows: 

d -§- = - E eg k + E v m ^ c+EEE"« ^ «r ^ - 

m— 1 m— 1 r— 1 p—1 n— 1 

4.37 

3 3 3 3 3 3 3 v ' 

E E E acm u mrs v r jf + E E E E 3— ^ mrs ^ ^ 

m— 1 r— 1 s—1 m— 1 r— 1 s—1 p—1 

Now let's apply the operator V r to the equality (4.35). As a result we obtain 

3 3 3 

w r jf = -w r ef + E E E "r-y* v ^ 7 9 ™ Ri + 

r=l7=1 " =1 (4.38) 

3 3 3 

+ EEE^ ^ g Tn V r Rt 

r — 1 7—1 n— 1 

Taking into account (4.38) and (4.36), we can perform the following calculations: 

3 3 3 3 3 3 

E E E 3am W mra V r = " Y, E E 3am ^ V r 0? + 



m— 1 r— 1 s—1 m— 1 r— 1 s—1 

3 3 3 3 3 3 



E E E E E E 9am 0J T1S W Q™ flg + 

rn— 1 r— 1 s—1 r — 1 7— 1 n— 1 
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333333 333 

+ E EEEEE2- ^ mrs ^ « 7 g Tn v r ji£ = - EEEff- x 

r— 1 s— 1 r— 1 7—1 71—1 m— 1 r— 1 s— 1 

3 3 3 3 3 

x v r 0? + E E E E E a am (c <5; - s™ s r T ) w 9 ™ r£ + 

m—l r— 1 t — 1 7—1 n— 1 
3 3 3 3 3 

+ E E E E E (c - s m s) « 7 <? rn v p j# = 

m— 1 r— 1 r — 1 7— 1 n— 1 
3 3 3 3 3 3 

= - E E E 9am " mrs v r ef + ^ - E E 5™ ^ + 

m=l r— 1 s— 1 r— 1 r— 1 n=l 

3 3 3 



E^' V r<"EE B » 5™ V rJ Rf. 



r— 1 r— 1 n— 1 

The result of these calculations is used in order to transform the equation (4.37): 

d -§- = - E ^ *c + E c+EEE"« ^ 3™ fl£ - 

m— 1 m—l r— 1 p— 1 n— 1 

3 3 3 

- iy r v r R a + v r V r R a ) + E 5™ ( V ^« + v « V r i^ ) + (4.39) 

r— 1 r— 1 n— 1 

+ E E E - mrs + E E E E - mrs ^ */■ 

m—l r— 1 s— 1 m—l r— 1 s— 1 p— 1 

The following identity is known as the zero divergency condition (see (4.8) in [3]): 

£ £ <r v r i^ + E E E ^ 3™ ^ = °- ( 4 - 4 °) 

r— 1 r— 1 p— 1 r— 1 n— 1 

Applying (4.40) to (4.39), we find that two terms in (4.39) are canceled: 

anfi 3 3 3 

= y: v m ^ - x ^ Ri~Y- r + 

m—l r— 1 r— 1 

3 3 3 3 3 

+ E E 9 rn v r v a ^ + EEEff- w ' nrs v ^ + ( 4 - 41 ) 

r— 1 n— 1 m — l r— 1 s— 1 

3 3 3 3 3 

+ E E E E 9am CO mrS Z% Of i?™. 
m—l r— 1 s — 1 p— 1 m—l 

Having completed the transformations of the differential equation (4.33), now 
we return to the formula (4.32). Substituting (4.41) into (4.32), we obtain 

3 / 3 \ 3 3 

d t e + X v k dkS — T I d t s + E ^' 9fcS + E E ^ af3 W x 

fe=l V fe=l / a=l/3=l 
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/ 3 \ 3 3 



V m=l ) a=l (3=1 

( 3 3 3 3 

E v ^ m < - E v ^ r R * + E E ff" 1 v ^ ^ + 
m— 1 r—1 r—1 n— 1 

3 3 3 3 3 3 3 

E EE^™ ^ mrs v^f + E E E E ^ ^ z% e? 



m— 1 r—1 s—1 m— 1 r—1 s—1 p—1 

3 \ 3 3 

- E C R rn + E E ^ <° ^ V «^' 

m=l / a=l/3=l 

Let's perform some rearrangement of terms in the above equality. We write it as 
3 / 3 \ 3 3 3 

d t e + Y / v k d k e = T Us + E^s +EEE T x 

fc=i V fc=i / 1=1 j=i fc=i 



(4.42) 

V fe (/'*•' %) , ^ ^ V, »j S« % 



p * — — ' p A — — ' p 

r i=i j=i r i=i j=i r 

For this purpose we introduce some auxiliary notations. We denote 

3 

% = E^/- ( 4 - 43 ) 

k=l 

It is easy to see that the formula (4.42) expresses the standard index lowering 
procedure. The other notations are less trivial. We denote 

g g g „ „0i n i ,mkj 

pkii - E E E t v »- ( 4 - 44 ) 

a=l (3=1 m=l 

The quantities (4.43) define a tensor field of the type (3, 0). Then we introduce the 
symmetric tensor field er of the type (2, 0): 

a H = p (V«[2]e + V 3 '*[2]e) + 

3 3 

+ E E P W £ R ° 9 ai + Vj[4]£ Ri g°*) - (4.45) 

a=l/3=l 

3 3 

-EE/ 3 W £ ^ + *° 2 ^ 

a=l ,3=1 

And finally, we introduce another tensor field of the type (2, 0) with the components 

3 3 

6« = p (V*'[2]e + V ji [2]e) + E E P V P^ £ R * 9™ ~ 

" = 1/3=1 

3 3 3 3 3 3 ( 446 ) 

- E E E E E P V e^ £ 9*m U™* Z% <r + j> V, /-''. 

en— 1 (3—1 m— 1 r—1 p—1 A:— 1 
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Note that (4.45) is not simply a notation. The same symbol a is used in (4.25) and 
in (4.3), where it represents the components of the regular stress tensor. 

Theorem 4.4. In the case of frozen dislocations the components of the regular 
stress tensor <j are given by the formula (4.45). 

Proof. Let's remember that the case of frozen dislocations is that very case where 
the dislocation lines move together with the medium like water-plants frozen into 
the ice (see [2]). In this case we have purely elastic response of the medium: 

0/ = (4.47) 
(see theorem 2.2 in [3]). From (4.43) and (4.47) we derive 

% = 0. (4.48) 
Let's substitute (4.48) into (4.42) and replace a 1 ^ by er u in this formula: 

d t e + J2v k d k e = T Us + £ v k d k s) + £ £ ^-M. (4.49) 

k=l \ k=l ) i=l 3=1 P 

Then we substitute (4.49) into (4.25). Taking into account (4.3), (4.4), and the 
first balance equation (4.1), from we (4.25) derive 



k=l \ i=l / i=l j=l 

1 12 12 12 J2 — — — [ 12J2 VtT — VkT 

i=l k=l j=l 9=1 ^ i=l k=l ^ 

where Vij are determined by the following formula: 



(4.50) 



(4.51) 



The equation (4.50) has a quite transparent interpretation. This is the entropy 
balance equation. Each term in its right hand side corresponds to some definite 
mechanism of entropy production. Last two terms correspond to viscosity and heat 
conductivity phenomena respectively. In the case of frozen dislocations, i. e. in the 
case of purely elastic medium response, we have no additional entropy production 
mechanisms. Hence, the first term in the right hand side of (4.50) should vanish: 

a ij _ frU = o. (4.52) 

The equality (4.52) means that a*i coincides with cr'-' , i. e. (4.45) is a true expression 
for the components of the regular stress tensor in the case of frozen dislocations. 
The theorem is proved. □ 

Conjecture 4.1. The components of the regular stress tensor cr for a dislocated 
crystalline medium are always given by the formula (4.45). 
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Assuming that the conjecture 4.1 is valid, we substitute (4.42) into the equality 
(4.25). As a result from (4.25) we derive the entropy balance equation 

fc=l \ i=l i=li=l / ^gg^ 

~ £ Z ~ g °' Vtk ^ V3q 1 Z Z v ' r VkT 



i=l j=l i=l fe=l j=l g=l i=l fc=l 

where and u, g are determined by the formula (4.51). 

5. Some conclusions. 

The equation (4.53) is a basic equation for understanding the thermodynamics of 
plasticity in crystals. It is similar to the equation (10.13) in [5]. However, there are 
some visible differences. In the left hand side of (4.53) we have the additional term 
produced by the tensor P with the components (4.44). It describes the entropy 
carried by moving dislocations. As for the first term in the right hand side of 
(4.53), it is also different from such a term in [5] since ^ 6 U (compare (4.45) 
and (4.46) above). The entropy growth condition leads to the inequality 

i=i j=i 

It is similar to (10.14) in [5]. Note that in the present theory we have no restrictions 
like det G = dct G and det G = 1. We got rid of them by including the density p 
as an explicit argument in (2.8) and (2.12). As a result we have no restriction for 
the trace of the matrix 0?. The symmetry condition for 9ij is also absent. 
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